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GENERAL CONSIDERATIONS ON THE FLOW OF COMPRESSIBLE FLUIDS* 

By L. Frandtl 



preliminary Remark 

The Royal Academy of Italy has conferred upon me the 
honor of opening the theoretical division of this year's 
Volta Congress "by a lecture which shall serve as an intro- 
duction to the subject of the flow of compressible media. 
In view of the limited time available, I shall entirely 
omit any description of the historical development of this 
branch of the science and shall also pass over the very 
familiar phenomena of hydraulics that receive the usual 
one-dimensional treatment as, for example, the discharge 
from orifices, the Laval nozzle, etc. I shall confine my- 
self, rather, to a consideration of the most important 
properties of these flows, from a modern point of view, 
starting from the differential equations of compressible 
flow . 

I. INTRODUCTORY CONSIDERATIONS 



The problem of the motion of fluids which is already 
sufficiently involved even when considered as incompressi- 
ble, becomes still further complicated and more difficult 
when the property of compressibility is taken into account. 
In the majority of cases, therefore, when the compressi- 
bility is to be allowed for, we are forced to make simpli- 
fying assumptions in some other direction. Thus in our 
discussion we shall have to neglect viscosity and so as- 
sume our fluid to be frictionless and compressible. We 
shall further assume that the density of the compressible 
fluid depends on the pressure only and such inhomogenei- 
ties as, for example, the heat conducted to the fluid from 
the outside or arising from inner combustion, are excluded 



* H Allgemeine Ueberl egung en uber die Stromung zusammen- 
druckbarer Flus s igke i t en . 11 (paper on High Speed in 
Aviation, read at the 5th Volta Congress, at Rome, 
September 30 to October 5, 1935.) Reale Accademia 
d 1 Italia, Classe delle Scienze Fisiche, Matematiche 
e Natural i. 
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from our discussion. We 
the pressure p and the 
mined . 



assume that 
density p 



the relation "between 
is uniquely deter- 



For this f rict ionless , homogeneous, compressible flu- 
id, just as for the f r i ct io nl e s s , homogeneous, incompres- 
sible fluid, the law of Lagrange applies, namely, that a 
fluid without initial circulation continues to move with- 
out circulation. It may be observed that this law holds 
true for steady motion only. Where the velocities are 
higher than that of sound, the motion begins to be un- 
steady with the occurrence of irreversible compressions, 
leading to an increase in entropy. In this way homogenei- 
ty is in general destroyed and the Lagrangian law is no 
longer appl i cable • 

Since the state of rest is a special case of irrota- 
tional motion, any compressible fluid flow starting from 
rest, whether steady or unsteady, is an irrotational flow, 
and as such may therefore be represented by a velocity po- 
tential whose gradient is the velocity. 



V == grad $ 



(1) 



the components "being 



3$ 
3x 



v = -- 



By' 



w = 

dz 



(la) 



For such potential flows of homogeneous frictionless 
fluids the Bernoulli equation may be applied: 



at 2 



30 



+ P - U = f(tj 



(2) 



In this equation 
the force 



p - / 



denotes the pressure func- 



tion and U 

is the force considered U 
function of time. Another 
the equation of continuity 
constancy of the mass. We 
that a definite element of 
fluid continues to maintain 



P 

function; in the case where gravity 
= - gz; f(t) is an arbitrary 
equation to be considered is 
which is an expression of the 
may express it either by saying 
volume of the same particle of 
a constant mass or that the 



mass of an element of fixed volume 
that more flows in than flows out. 
view we arrive at the equation 

|| + div ( P V) 



changes with time so 
From either point of 



(3) 
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where it is possible to put div (pV) = p div V + V • 
grad p (where the dot denotes the scalar product). Equa- 
tions (2) and (3), together with (l), define our problem. 
Except where meteorological applications are considered, 
gravity generally plays a subordinate part in problems 
connected with compressible flow, so that it is always 
possible to neglect the term containing U. In such cases 
as are considered in acoustics where the velo c i t ie s g ar e 
small and rapidly changing, the quadratic terms (V /2 in 
equation (2) and V • grad p inequation (3)) may be neg- 
lected. By equation (l) and the relation div grad $ = A$ 
equation ( 3) may then be written 

|~ + pA$ = 0 (3a) 



We now differentiate 
and write 



equation (2) with respect to time 



3? 

at 



dp 

dp 



9P 
3t 



sp- 
at 



(2a) 



The dimensions of dp/dp are those of a velocity 
squared and it may therefore be set equal to c2 , where 
the velocity c is still a function of p; when the ve- 
locity is zero at infinity, f(t) may be set equal to a 
constant. The differentiated eouation (2) thus becomes 



6 2 $ 
9t 2 + 



p 9t 



(4) 



Prom the relations (3a) and (4) the equation is further 
reduced to 



aT 2 ' = c 



(5) 



which is the familiar differential equation of sound. A 
well-known solution of this differential equation is the 
plane sound wave, corresponding to the eqiiation 

$ - F ( x - ct) (6) 



and which travels with velocity c, the form of the wave 
being given by the quite arbitrary function F. Another 
solution is the spherical wave whose equation is 



J 
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$ = L f ( r - ct ) (7) 

. ' r • • 

where r denotes the distance from a center. The poten- 
tial for the cylindrical wave is less simple. It may "be 
written 

^ 2 f (h dl 
■■'•*/ (8) 
i x /(ct - If - F 

In all these examples c clearly signifies the veloc- 
ity with which the waves travel outward. It is customary 
to call this velocity, which we have assumed above as a 
function of the density, the velocity of sound. The fact 
that in the above equations we have treated it as a con- 
stant clearly shows that the solutions are only approxi- 
mations for very small sound amplitudes where the density 
nowhere differs much from the initial density p Q . For 

finite amplitudes it is necessary to go hack to the exact 
equations (2) and (3) and we thus arrive at the result 
first obtained "by Riemann * that the regions of condensa- 
tion move forward at a greater velocity than the regions 
of rarefaction and therefore the wave front "becomes steep- 
er and steeper until an unstable wave is formed which is 
known as a 11 shock wave. 11 Since in the condensed region 
the gas moves forward, and moves "backward in the rarefied 
region, while the impact wave moves with velocities that 
do not differ much from the sound velocity, we have the 
result that the relative velocity between the shock wave 
and the condensed region is smaller, and that between the 
impact wave and the rarefied region larger, than the sound 
velo city, 

II. STEADY POTENTIAL FLOW 



The above qualitative discussion will be sufficient 
for the case of the unsteady sound motions with small flow 
paths. I shall now turn to the main subject of my lecture, 
namely, steady potential flow. If we again neglect grav- 
ity equations (2) and (3) assume the form 

v 2 

~ + P (p) = const (2b) 

and 

div V + - V • crad p = 0 (3b) 



*See reference 1. 
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We may now set 

1 1 dp 1 

- grad p = - grad p = pr grad P 

and find grad P from (2^) . T7e thus obtain 

2 

div V - -g- V • grad ~- = 0 (9) 

Expressed in rectangular components u , v, w, we may write 
V • grad V /2 as follows: 

V . g rad £ = T • 4 T • T = »• £ ♦ ,» | ♦ | 

{Btx , 8v\ /9v 3w\ , /8w 3u\ 

+ uv W + 9^ + Uz + 3y /' + WU + dzV 

Equation (9) expressed in cartesian coordinates thus as- 
sumes the form 

3u / u 2 \ , 8v / v 2 \ 3w / w 2 \ 
§i V " c*J + By [} ~ c 2 ,' + 3z \ " c 2 J 

2 ( 8u dwV, Sw\ 

" c 2 V* v 3y + vw 3x"-' + wu 3x7 = 0 

du dv 

(In this equation u-je was made of the relation g-— = g£ 

derived from the condition of irrotational motion, and the 
corresponding relations for the other components.) 

From the form of equation (9) or (9a) it is immedi- 
ately evident that in all cases where the velocity compo- 
nents arc all small in comparison witll the sound velocity, 
the relation reduces to 

liV 7 a 0 

o r 

du + & v + 9 w s q 
3x dy 3z 

When the resulting velocity is not van! shingly small com- 
pared to the velocity of sound, tut still is smaller than 
the latter, then the quantitative relations are changed 



(9a) 
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although the character of the motion still remains entire- 
ly similar to that of incompressible fluids. To see what 
equation (9) tells us when the flow velocities are of the 
order of magnitude of the sound velocity, we may so ar- 
range the system of coordinates that the X axis coincides 
with the direction of flow. Then only the velocity compo- 
nent u, in the neighborhood of the origin, will "be of 
the order of magnitude' of c, while the velocity compo- 
nents v and w will "be small compared to c. Equation 
(9a) then reduces to 

|s ( % m 4) + |v + |w (10) 

3x \ c 2 / 3y dz ' 

If we now assume that u increases in the direction 

of the flow corresponding to a fall in pressure in that 

direction, then 3u/dx is positive. From equation (10), 

dv Qw o p 

therefore, ~ + is negative when u < c and posi- 

dy 6 z ^ 

tive when u 2 > c 2 and becomes zero for u 2 = c 2 . Ex- 
pressed in words, this means that the flow converges with 
increasing velocity when the velocity is less than that 
of sound but diverges when the velocity is greater than 
that of sound, and moreover, that when the velocity pass- 
es through the sound velocity the streamlines are paral- 
lel. This conclusion fully agrees with what is found for 
the flow through a Laval nozzle from elementary considera- 
tions. In that case, too, it is found that when the pres- 
sure falls throughout the length of the nozzle the veloc- 
ity in the converging part of the nozzle is smaller than 
the sound velocity and in the diverging part is greater. 
The sound velocity is exceeded just where the cross sec- 
tion of the nozzle is a minimum. 

Equation (10), however, tells us more. If we intro- 
duce the relation given by equation (l), we obtain 

d 2 $ / u 2 \ d 2 $ 9 2 $ „ , s 

If we consider the immediate neighborhood of some def- 
inite point, it is sufficient to set u equal to its mean 
value in the region under consideration. Using this sim- 
plification, there is obtained a very well-known partial 
differential equation. This is of the elliptic type for 
the case where the expression 1 - u 2 /c 2 is positive; 
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that is, it is entirely related to the differential equal 
for incompressible flow. When 1 - u 2 /c 2 is negative, 
however, which is true when the velocity u is greater 
than the sound velocity, then the elliptic type goes over 
into the hyperbolic type. It is known, however, that the 
solutions of the elliptic type are regular within the re- 
gion, whereas those of the hyperbolic type also admit of 
discontinuous solutions within the region, running through 
the range of the so-called characteristics of the differ- 
ential equation. What is typical of "both cases may he 
seen from the following consideration, which starts out 
from the well-known singularity in a flow source of an in- 
compressible fluid 



Let us inquire what is the form of the potential of 
such a source when the flow from the source has a constant 
velocity of the order of magnitude of the velocity of sound 
superimposed upon it. The problem may be simplified by 
assuming that the velocities from the source are small, 
thus making equation (10) or (10a) applicable. To find 
the singularity for this case we may apply a relative sys- 
tem in which the undisturbed medium is at rest and the 
source moves with the constant velocity ~u Q . We may first 
assume an 11 explosion wave" (knallwelle) of spherical shape 
expanding in all directions, the potential of which is ob- 
tained by assuming in equation (7) a function f, which 
is different from zero only within a very short interval 
and vanishes outside this interval. This is the case when 
a small volume suddenly begins to increase and then main- 
tains its new magnitude. It may then be assumed that a 
continuous series of such short expansions proceeds in 
such a manner that the center of the expansion travels 
forward with the velocity -u 0 . Since we are here consid- 
ering sound wave expansions obeying the differential equa- 
tion (5), the potential for the whole process may, on ac- 
count of the linearity of the differential equation, be 
built up by the superposition of the potentials of the in- 
dividual waves. When the expansion has continued long 
enough and the center of the expansion is momentarily at 
the origin of coordinates, the computation leads to the 
f o 1 lowing f ormula : 




A 



r 




A 



(11) 
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When u 0 is smaller than the sound velocity c, then 
the surfaces of constant potential instead of "being spher- 
ical as is the case for incompressible fluids, are flat- 
tened ellipsoids of revolution. The more closely the ve- 
locity u 0 ' approaches the velocity of sound the more 
strongly flattened the ellipsoids "become. When u 0 ex- 
ceeds the velocity of sound, however, then the solution 
is different from zero only within a cone of angle a, 
which is determined by 



tan a = 




or sin a = ± ~- (12) 

u o 

The same result is also obtained when the momentary 
position of the individual waves is investigated (fig. l) • 
The angle 'a is known as the Kach angle. (it may here he 
remarked that when u Q < c* the waves fill space in all 
directions, whereas for the case u 0 > c, they fill only 
the cone of fig. 1.) For u 0 > c, the constant potential 
surfaces are hyperboloids of two sheets having the given 
cone as asymptotic surface. Only one of these sheets has 
physical reality, as is easily seen from our description 
of the formation of the source potential from the explo- 
sion waves. Where still another potential is to be built 
up by the superposition of such sources, it is to be ob- 
served that everywhere outside the given cone, the poten- 
tial of the individual sources is to be set equal to zero 
and takes on the values of formula' (ll) only within the 
cone. 

III. PLOWS WITH 11 SUBSONIC 11 VELOCITIES - LINEAR THEORY 



Formula (ll) enables us to obtain a general relation 
for the solution of the differential equation (10a) and 
this we shall discuss further on. For "subsonic" veloci- 
ties it is possible to relate every solution of this dif- 
ferential equation to a solution of the differential equa- 
tion for constant volume potential flow 

2 2 2 

9 $ 9 <£ 9 $ 

W + a? + ap a 0 
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and assign the same potentials to the corresponding points 
of each space (reference 2). The assumption must naturally 
he made, as for differential equation (10a), that the ve- 
locities derived from the potential are small compared 
with the "basic velocity u G . The application cannot there- 
fore "be made to flows in which a stagnation point occurs, 
since in this case the deviation from u Q is just as 
large as u 0 itself. It is permi s si ble , however, to ap- 
ply the .equation to flow about very narrow shapes having a 
sharp entrance edge. Given the potential flow $ 0 for an 
incompressible fluid, the question next arises as to what 
is the form of the contour that corresponds "by the above 
rules to the potential of the compressible fluid. The 

slope ~ = ^ may, to a sufficient approximation, "be put 

equal to 

v _ JL_ 30 
u u 0 9y 

Co rre spondingly , 



dri 1 30 8 _ d$ 
' - — ; with -- 

<H u 0 dr\ 9y 




we have — - -j / 1 - -j— 



dx d£ 



c 



We then obtain the result that the entire contour must "be 
made thinner than that corresponding to the incompressi- 
ble fluid with equ al potential values approximately in the 
/ u n 2 

ratio* J 1 - and similarly the angle of attack must 

*It may "be pointed out that the points on the contour do 
not correspond to the above transformation equations for 
the potential, which require an increase in the y ordi- 
nates for points of equal potential. The above conclu- 
sions are nevertheless applicable since the difference in 
the flow direction for points lying near each other differ 
only by a second degree order from the result given above 
on account of the slenderness of the contour under consid- 
eration. 
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"be made smaller in the same ratio.* 

We must still consider the pressure distribution on 
the surface. This is determined according to the Ber- 

3u 

noulli equation by the tern pu g-- for which we may put 
pu 0 ~f- . Since we had set % = x and $ was to have 

the same value at corresponding points, this magnitude re- 
mains constant in the transformation, and the same holds 
true for the pressures themselves. The tendency toward 
flow separation may therefore also he expected to "be the 
same for both cases. The conclusion follows that in or- 
der to avoid separation of flow it is necessary to make 
the profiles of the airfoils, etc, , flatter and the angles 
of attack must he made correspondingly smaller, as the ve- 
locity u Q approaches the velocity of sound. The maxi- 
mum lifts attainable according to this approximate theory 
will he the same as those that may he expected for the in- 
compressible flow. Actually, it is foundihat on approach- 
ing the velocity of sound, the relations are considerably 
less favorable than indicated by this approximate theory. 
The chief reason for this is the fact that the superim- 
posed velocities are not actually small compared to the 
basic velocity u 0 and consequently it is possible for 
points on the suction side to exist at which the velocity 
of sound is either reached or exceeded, so that considera- 
ble deviations from our computations are to be expected. 

It may still be asked how the fact is to be explained 
that as the velocity of sound is approached, the same lift 
is obtained with a less cambered profile and at a smaller 
angle of attack. The explanation lies in the fact that 

ill 

It is also possible to coordinate the points in such a 
manner that the potentials 0 in the xyz space will be a 
multiple of the potentials at the corresponding points of 

the lr\t space. If the multiple chosen is l/ / 1 - ~& — 

then ~ = that is, the profiles and angles of attack 

dx 3| 

will now agree. The differences in pressure (see below) 

will now be raised in the ratio if 1 1 - ~- and the 

v c 2 

tendency to separation will therefore be increased. 
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according to the transformation given above, the potential 
field and therefore also the velocity field in the direc- 
tion at right angles to the "basic velocity extends further 



extended region a total impulse of the same magnitude as 
would "be the case for incompressible flow. 

For "basic velocities above that of sound this analogy 
cannot naturally "be applied, since the transformation for- 
mula (13) would in this case give an imaginary result. 
Here other methods must he used. (See sections 5 and 6.) 



IV. FLOWS WITH SUBSONIC VELOCITIES - HIGHER APPROXIMATIONS 



If it is desired to obtain a "better agreement with 
fact than is afforded "by the linear theory, it is possible, 
starting either from the theory of potential flow for in- 
compressible fluids or from the solutions of the differen- 
tial equation (10a), by a step-by-step process to obtain 
closer approximations resulting from the application of 
the exact equation (9a). Computations of the first kind 
were carried out by many investigators (references 3 and 4) 
but, to the best of my knowledge, no computations of the 
second kind.* In practice these computations are rather 
laborious and do not agree very well with each other. It 
is therefore a noteworthy fact that it is also possible to 
find good approximations on the basis of electrical anal- 
ogy, whereby a good approximate solution may be obtained - 
only for the case, however, where the velocity of sound is 
not attained at any point of the region considered. Since 
this method has been worked out by Professor G. I. Taylor 
(reference 7), who will give a more detailed report on it 
at this session, I shall not touch on this subject any 
further . 

Through three short notes by Riaboi2chinsky (reference 
8) and Demtchenko (references 9 and 10), in the Comptes 



which are smaller in the ratio 



in the ratio 




so the vertical velocities, 



produce in this 



*&• Braun (reference 6) has obtained solutions in the 
neighborhood of the sound velocity by applying a variation 
principle of Bateman (reference 5), 
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Rendus of the Paris Academy, 1932, and old work "by 
Tchaplygin (reference 11) written in the Russian language 
in 1904 "became known, wherein it was shown that for the 
case of two-dimensional flow, the problem may, "by trans- 
formation to new coordinates (rectangular or polar), he 
presented in such a form that a potential flow "between two 
plates at a predetermined variable distance apart, may he 
represented by an electric flow in an electrolyte of vari- 
able depth. One coordinate x will be a function of the 
ratio of the local velocity to the maximum (which may nat- 
urally also be written as a function of the ratio of the 
density to the maximum density); the other coordinate P 
is the direction angle of the velocity of the flow, the 
variable distance between the two plates being a function 
of the first coordinate only. Since it is possible to 
compute mathematically the solutions for such regions 
whose boundaries are determined by the lines X - constant 
and 8 = constant, it becomes possible to solve those 
problems of the Helmhol t z-:'Iirchhof f type where flat walls 
(direction 6 = constant) and free boundaries (constant 
velocity, therefore X ~ constant) are considered for the 
case where compressibility is taken into account. For one 
particular condition which, though not occurring in na- 
ture, may be applied approximately when the difference in 
density is moderate, the problem nay even be formulated 
as a usual type of potential flow in the coordinates X 
and 6. Bemtclxenlco (reference 10) showed among other 
things that for the Eirchhoff flow against a flat plate, 
when the velocity of the flow is half that of sound, the 
resistance is about 7 percent higher than the Kirchhoff 
resistance. 

V. PLOW WITH SUFERSOKIC VELOCITY - LINEAR THEORY - ELOW 



I shall now consider the flows for which the velocity 
is greater than the velocity of sound (the so-called "su- 
personic velocity"). We may here again start out from 
formula (10) or (10a). If a two-dimensional flow is con- 
sidered, then equation (lOa) now reads as follows: 



AROUND A CORNER 




(10b) 



The general solution of this linear differential 
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equation may immediately "be set down. It is 



0 = Pi (y - x tan a) + F 2 (y + x tan a) 



(14) 



where again 



tan a = l/ 



/ 




(15) 



so that 



c 



sin a 



o 



According to the reasoning that led from equation (9a) 
to equation (10), equation (14) is the solution for the 
neighborhood of a point, when the X axis is chosen to lie 
in the direction of the flow. The functions F x and 3? a 

are subject to the limitation only that the magnitudes ap- 
pearing in the differential equation (10b) exist; otherwise 
they may be taken quite arbitrarily. When the functions 
F x and F 2 are taken to represent a wave form, then equa- 
tion (14) represents the superposition of two wave trains, 
crossing each other at the angle ±a with respect to the 
mean direction of flow. For a more general solution, the 
wave trains are such that for each position represented by 
equation (15), u Q corresponds to the local mean veloc- 
ity at that point. Since such waves may be photographed 
(see fig. 2), it is possible by measuring the angle 2a 
between the two crossing wave trains to obtain the ratio 
c/u 0 . (It should be noted that c is not the constant 
sound velocity of the gas at rest but the variable sound 
velocity of the adiabat ically cooled gas.) The direction 
of the flow bisects the angle. 

Let us now consider the following special case of 
flow. In equation (14) let F 2 = 0 , F x = 0 for positive 
values of the argument y - x tan a; for negative values 
let F x = A (y - x tan a). In the neighborhood of zero 
it is possible to pass to the limit zero by means of a 
transition arc, or the radius of curvature of the transi- 
tion arc may be used. The above formula corresponds to 
the limiting value. To obtain the total potential, the 
potential of the straight streamlines u 0 • x must be 
added. In the region x > y cot a the components of the 
velocity then become 
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X tan a 



a y ~ X 



(16) 



In our example, therefore, the flow is parallel in 
front of the boundary line y = tan a and is in the direc- 
tion of the X axis. Behind the "boundary line it is also 
rectilinear hut makes an angle |3 to the X axis given by 
(see figs. 3 and 4) 

r> v X X / N 

tan 8 = - = r— r (17) 

u u n — A tan cc u n 



If the density and pressure upstream from the "boundary 
y = x tan a are denoted by p 0 and p 0 , respectively, 

and the pressure "below the "boundary "by p, then approxi- 
mately 

P - Po = - Po u o ' (u - u o> 

(18) 

= p 0 u Q X tan a 



Along the boundary y = x tan a, therefore, there is a 
pressure jump (which may he converted into' a steady pres- 
sure rise "by introducing the transition curves mentioned 
above in the function Pi). The fluid particles are accel- 
erated in the direction of the pressure rise, that is, nor- 
mal to the boundary line y = y tan a. This result may 
also be obtained directly from equation (16), from which 
the direction of the velocity vector is obtained, as fol- 
lows: 

tan a' = -J. = - 1 

u - u 0 tan a 

Positive values of X and (3 correspond to a conden- 
sation; negative values to a rarefaction. It should be 
stressed, however, that the above computation is derived 
from the approximated linear differential equation (10a) 
and is therefore valid for small pressure differences and 
small angles of deviation only. 

What may be expected for larger deviations may clearly 
be seen if we take two pressure jumps of the same sign at a 
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short distance from each other and see what we get from 
equations (15) to (18). The result is obtained that for 
the case of a pressure rise the two pressure jumps con- 
verge; for a pressure drop they diverge. Not only is the 
direction of flow in the second field rotated "by the an- 
gle p with respect to the first, hut the velocity "behind 
the first pressure wave is smaller and therefore the Mach 
angle oc 2 is larger than CCi. For a rarefaction wave, on 

the contrary, the velocity hehind the first wave is larger 
and therefore the l.Iach angle smaller. The two wave fronts 
converge or diverge accordingly "by the amount (3 4- CC2 - CCi. 

The transition from the type of flow just considered 
to that of a continuously curved flow may now easily he 
effected hy replacing the first streamline (given wall 
al-ong which the flow takes place) hy a polygon of very many 
sides. From each edge of the polygon, waves start out ex- 
panding and diverging in case of a convex wall and con- 
tracting and converging in case of a concave wall. In 
this latter case it is also possible that the waves meet 
each other completely. The compression may then assume an 
unsteady finite value (the so-called "compression shock"). 
This behavior of compression and condensation waves start- 
ing- out from a curved wall is entirely analogous to the re- 
sult already mentioned obtained by Riemann for plane waves 
of finite amplitude. 

We see, moreover, from the result of a transition to 
a boundary of radius of curvature zero, that for a wall 
which forms a convex side inclined at an angle -|3 (see 
fig. 0), the flow remains unchanged up to a surface form- 
ing an angle oci with the direction of flow, then expands 
within an angle -p + CCi - oc 2 maintaining the pressure 

and velocity constant along each ray; then in the direc- 
tion -|3 again passes over into parallel flow with con- 
stant velocity. The quantitative relations for a flow of 

this type for a gas that obeys the law P = P 0 • (p/p 0 ) 

were given in the Gottingen dissertation of Th. Meyer in 
1908.* In case of a concave edge there is a compression 
shock which lies between the angles oci and (3 + cc 2 . In 

front of and behind the condensation shock the flow shows 
constant velocity and constant pressure. 

*See reference 12 # 
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VI. APPROXIMATE METHOD FOR GENERAL TWO-DIMENSIONAL 

SUPERSONIC PLOW 



I shall now turn to the case where the flow is af- 
fected not only by a single wall but also by an oppositely 
lying wall. In this case the waves pass through each oth- 
er from opposite sides and also deviate from their origi- 
nal direction. The relations are most clearly "brought out 
by the following procedure. 

In the same manner as, for the purpose of certain 
approximate calculations, a curve is replaced "by a series 
of straight- 1 ine steps or a polygon, it is possible to re- 
place the continuous deviation of the velocity direction 
in a supersonic flow by a series of sudden deviations of 
the kind shown in figures 3 and 4, If the angle (3 of 
these deviations is chosen to be of the same size in all 
cases, for example 2°, then only such directions occur 
7; hich differ by an integral multiple of that is, in 

our example, of 2°. If we start from a definite super- 
sonic velocity, then on the basis of the Meyer formula, 
only certain discrete values of the velocity may occur. 
In order to obtain a better idea of these velocities and 
directions, the method suggests itself of drawing all the 
possible velocity vectors of the system at each point. If 
in this figure we trace those conditions corresponding to 
the case we have just considered, namely-, where waves are 
assumed to start out from one wall, then all the points 
in the velocity picture corresponding to that state lie on 
a single curve (the thick curve of fig. 7). The entire 
system of velocity vectors for the case where waves travel 
into the region under consideration from both sides is 
then evidently obtained by drawing curves of the same sort 
tiirough all points lying 4° apart on a circle w = constant 

(w = v /u 2 + v 2 ). All the intersection points so obtained 
then represent the end points of the velocity vectors possi- 
ble in this system. 

According to a relation shown in a previous section, 
the direction of the vector difference of the velocity, in 
front of and behind a wave, is perpendicular to the wave 
front. This relation enables us to construct graphically 
the entire flow picture for any flow when the magnitude 
and direction of the velocity at the entrance section, and 
cither the pressure or direction of flow at the side edges, 
are given beforehand. Por the permanent gases the curves 
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of figure 7 are epicycloids and it is possible to construct 
the curves easily on the drafting "board "by using a rotata- 
ble ellipse to draw the directions normal to the curve tan- 
gents, as Busemann ha-s shown. 

In order to maintain constant the flow direction along 
the wall, for the case of a straight wall, or maintain the 
velocity constant in the case where a constant pressure is 
given, in accordance with Bernoulli's law, it is necessary 
to reflect the expansion and compression waves that reach 
the "boundary. As may easily "be seen from figure 7, at a 
straight wall condensations are reflected as condensations 
and rarefactions as rarefactions; at a free jet "boundary, 
however, condensations are reflected as rarefactions and 
rarefactions as condensations. For curved walls, the 
walls are replaced "by a polygon of angle (3 and waves 
start out from each corner. 

In order to obtain a clear construction it is advisa- 
ble to number each epicycloid, one set with numbers increas- 
ing in clockwise order, the other set decreasing in clock- 
wise order. It may then easily be seen that the difference 
of these sets of numbers is a measure of the angle of the 
flow direction, while the sum is constant on a radius; in 
other words a function of the velocity. It may also be re- 
marked in this connection that, as in the treatment of the 
subsonic velocity by Tschapligin, the angle of the flow di- 
rection and a function of the velocity again appear as de- 
ciding factors. 

It is possible by a contact transformation to trans- 
form equation (9a) so that u, v, and w become the inde- 
pendent variables; and in the case of two-dimensional flow 
in the plane of u and v, it is possible to pass over to 
polar coordinates and so introduce the radius and the an- 
gle as new independent variables. As the form of equation 
(9a) already shows, the contact transformation offers the 
advantage that the differential equation becomes linear in 
the dependent variables, since the expressions 1 - u 2 /c 2 , 
etc., are now functions of the coordinates. Ey a suitable 
stretching of the radii - that is, by introducing a func- 
tion of the radius instead of the radius - it is possible 
to simplify still further the differential equation as was 
done by Tschapligin. When the relations between the waves 
in the field of flow and the corresponding velocity field 
are somewhat more closely analyzed, it is found that to 
each field of flow limited by two waves in the one system 
and two waves of the other system there corresponds a point 
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in the velocity field (since according to the assumption 
made the velocity is to "be constant in each such field). 
Conversely, to a point of intersection of two waves in the 
field of flow there corresponds a quadrilateral in the ve- 
locity field, since at the point of intersection four 
fields of different velocities meet. Both of these fig- 
ures are therefore the reciprocals of each other, just 
like a truss diagram and its reciprocal force diagram. 
The similarity goes even further. Just as a truss togethBr 
with its force plane may be represented "by the Airy stress 
surface, so in this case there exists a corresponding 
function, namely, the function ux + vy - $> obtained from 
the potential 0 "by the contact t ransf ormat ion • 

The graphic method here described, to which there nat- 
urally corresponds an analytical method derived on the same 
basis, offers a very convenient means of discussing all two- 
dimensional supersonic flows. Applications of this method 
have been made to the discharge from nozzles, the flow 
about streamlined bodies, airfoils, etc. See A. Busemann 
(reference 13) or the examples given in his article in the 
Kandbuch der Experimentalphysik (reference 14). For rota- 
tional symmetric flows there are no such simple methods. 
The analogous cases to the flow around a corner, the flow 
at the tip of a cone of finite angle, have been solved 
graphically by Busemann (refarence 15) and analytically by 
G. I. Taylor and J. W. Ma ceo 11 (reference 16). Solutions 
for the general case of supersonic flow about a body of 
rotation have been given at least by approximate linear- 
ized methods by Th. v. Karman and N, B. Moore (reference 
17) , which methods are chiefly applicable to very slender 
bodie s . 

VII. PASSAGE THROUGH THE VELOCITY OE SOUND 



The methods described in sections 3 and 4 for the 
flow with subsonic velocity approach one another only 
provided the velocity of sound is nowhere reached. The 
methods discussed in sections 5 and 6 for supersonic ve- 
locities are similarly applicable only when the velocity 
is everywhere greater than the sound velocity. Neither 
of these two methods therefore provides any information 
about the manner in which passage through the velocity of 
sound takes place. There are, however, solutions for 
such cases that have been obtained some time ago. The 
first one, given by T. Meyer (reference 12), was obtained 
by a power series development in the coordinates x and y 
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of a point lying on the curve u = c . It turned out that 
the flow potential could "be developed at such a point and 
tests showed good agreement with the computations. An- 
other example was given "by G. I. Taylor (reference 18), 
This is the curved line flow in a usual potential vortex, 
where the velocities are exactly as large as the velocities 
in the potential vortex of an incompressible fluid. Both 
of these examples, the one with a straight mean streamline 
and an increasing velocity in the direction of flow, the 
other with curved streamlines and a constant velocity along 
each streamline, allow us to conclude that in going from 
subsonic to supersonic velocities a continuous passage 
through the sound velocity may he expected. It is to he 
noted, however, that in the supersonic velocity region it 
is very easy for conditions to he set up such that a se- 
ries of infinitesimal condensation waves combine to form a 
condensation shock of finite magnitude. In the passage 
from large velocities to smaller such condensation shocks 
occur almost regularly. In this connection it should also 
he mentioned that, in contradistinction to potential flows 
with subsonic velocities, which are symmetrical forward and 
aft about bodies that are themselves symmetrical at both 
ends, a potential flow with supersonic velocities is proba- 
bly never symmetrical because the condensation and rarefac- 
tion waves starting out from the walls of the body always 
travel obliquely toward the rear and never forward. Per- 
haps it is still possible by taking this lack of symmetry 
into account, to solve analytically more of those cases of 
flow where passage from subsonic to supersonic velocities 
occur, especially when the supersonic region is limited in 
extent, as may be expected, for example, in the flow about 
cylinders, etc. 

From the analytical standpoint, for steady friction- 
less flow, it is possible by reversing the velocities of 
one solution to obtain another solution. In applying this 
method to a symmetrical body or channel with nonsymmetrical 
flow forward and aft, we thus obtain by reflection a sec- 
ond solution different from the first and with the same 
flow direction. The passage from the single valued sym- 
metrical solution for the subsonic velocities to the two- 
valued solution indicates the existence of a special kind 
of branching position. This may perhaps be the reason for 
the fact that the analytical methods that have been used up 
to the present have ceased to converge. By the application 
of the direct methods of the calculus of variations, G. 
Braun (reference 6) was able to obtain solutions for the 
branching positions, which solutions, however, should be 
applicable to infinitely small deviations only. 
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Here I should like to say a little more about the com- 
pression shocks. The points in front of and "behind the 
shock do not lie on an adia"batic, cut rather there is an 
increase in entropy due to the irreversible process. The 
magnitude that remains constant throughout the condensa- 
tion shock is the total energy = kinetic energy + heat 
content; per unit mass = w®/2 + I (the heat content I = 
U + pV is also called the enthalpy). In the general case 
the direction of flow forms an angle with the normal to 
the impact surface. These relations assume their simplest 
form for the two-dimensional compression shock. If the 
normal to the impact plane is taken as the X axis and the 
direction of flow is in the XY plane, the following equa- 
tions re sul t : 



Pi u i ~ P2 u 2 (continuity) 

2 s / 
$ + q 1 u 1 - p 2 + p 2 u 2 ° (impulse in the X 

direction) 

▼ l = v 2 . (impulse in the Y 

direct ion) 

| (u 1 2 +v 1 2 ) + l a = I (u 2 ^+v 2 2 ) + I 2 .. (energy) 



(The index 1 for the condition "before the impact, index 2 
for the condition after the impact.) 

The computations for this problem, using the equa- 
tions of state for ideal gases as a basis, were carried out 
by Th. Meyer (reference 12) after Ilugoniot (reference 19), 
and independently Stodola (reference 20) had previoxisly 
explained the behavior of the normal compression shock. 
The relations are clearly brought out in a diagram given 
by Busemann (reference 21). In the velo c it y- f ie Id picture 
he draws for each given initial velocity, the 11 impact po- 
lar , 11 namely, the geometrical locus of all the velocity 
vectors for the state after the impact and thus obtains 
curves like those of figure 8, where (3 denotes as before 
the deviation angle of the flow, y the angle between the 
impact plane and the direction of flow. The smallest val- 
ue of 7 is obtained for shocks with small velocity differ- 
ence and naturally agrees with the Kach angle a. The 
largest difference in velocity is obtained in normal im- 
pact. The velocity behind the normal impact is always 
smaller than that of sound. The increase in entropy prac- 
tically denotes a lowering in the pressure compared to 
that which would follow from the Bernoulli equation. This 
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difference in pressure is only insignificant, however, when 
the difference in velocity in front of and behind the shock 
is not too great. The process in any case involves a loss 
in energy. 



Before concluding we shall make a few observations on 
the application of the foregoing discussion to airfoils. 
As far as subsonic velocities are concerned, it may "be 
pointed out that in the neighborhood of the wing tip consid- 
erable supersonic velocities are set up which for high 
lift coefficients attain double the values of the flight 
velocity. It may therefore be expected that at speeds of 
170 m/ s (380 rniles/hr.) the velocity of sound may already 
be attained locally. This may explain the fact that in 
the region of 200 rn/ s (450 miles/hr.) there is already a 
notable decrease in the lift coefficient. 

As regards the induced velocity, it should be pointed 
out that it is possible to obtain the wing lift as well as 
the induced drag from the trailing vortices behind the 
wing. The velocities in these vortices are, however, in 
each case small compared to the velocity of sound, so that 
the usual laws for incompressible flow may be applied with- 
out objection for the computation of the lift and induced 
drag. In computing the lift distribution of a wing of 
given form, however, it is necessary to take into account 
the fact that due to the compressibility, the value of 



For the profile characteristics at supersonic veloci- 
ties, an approximate formula may be derived from the appli- 
cation of the considerations of section 5. If we consider 
a flat plate that is inclined by the angle £ to the di- 
rection of flow (fig. 9) (Ackeret reference 22), then on 
the pressure side the relations of figure 3 hold and there 
results therefore a pressure rise of an amount that is 
easily computed from formulas (l?) and (18). On the suc- 
tion side there is a corresponding lowering in pressure 
as required by the abo vo-ment ionod formulas for a negative 
angle. We thus obtain the lift coefficient 
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AFPLICATI01T TO AIRFOILS 



U(jq • 7-i 

is increased in the ratio 1 : / x - 

da v a a 

the approximate formula remains applicable). 




(as far as 




(19) 
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It is a noteworthy fact that the suction force, which 
exists in the case of subsonic velocities at the forward 
edge, is entirely missing in a flow of this -type and there- 
fore the resulting force is here not perpendicular to the 
direction of flow, , "but normal to the surface, so that even 
for a frictionless fluid the lift-drag ratio "becomes 

c c m fu 2 

It may "be seen that for a given value of c a this ratio 

is favorable just above the velocity of sound, hut with in- 
creasing velocity, becomes considerably less favorable, (At 
very large speeds the lift coefficients will likewise be 
very small.) These increased drags above the velocity of 
sound bear a connection with the waves that travel outward 
from the moving object. Buscmann (reference 23) was able 
to show that this wave energy was converted into heat part- 
ly near the wing and partly at a great distance away. 

As far as the induced drag is concerned, it follows 
different laws at speeds above the velocity of sound from 
those below the velocity of sound, Tor an unwarped rec- 
tangular airfoil the induced drag is all included in the 
triangular regions at each end that are limited by the 
Mach angle (cross-hatched areas in fig. 10). It is, more- 
over, not clearly separable from the wave resistance which 
according to the above is also proportional to c a 2 and 

therefore does not have the same independent significance 
as the case of subsonic velocities. According to a remark 
of Eusemann, the induced drag disappears entirely when the 
two triangular parts lying in the region of disturbance 
are removed. 



Translation by s. Eeiss, 
National Advisory Committee 
for Aeronautics. 
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